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1 Abstract 

2 At present, satellite remote sensing of coastal water quality and constituent 

3 concentration is subject to large errors as compared to the capability of satellite sensors in 

4 oceanic waters. In this study, field measurements collected on a series of cruises within 

5 U.S. southern Middle Atlantic Bight (SMAB) were applied to improve retrievals of 

6 satellite ocean color products in order to examine the factors that regulate the bio-optical 

7 properties within the continental shelf waters of the SMAB. The first objective was to 

8 develop improvements in satellite retrievals of absorption coefficients of phytoplankton 

9 ( a p h ), colored dissolved organic matter (CDOM) ( a g ), non-pigmented particles (ad), and 

10 non-pigmented particles plus CDOM (a</ g ), and chlorophyll a concentration ([Chla]). 

1 1 Several algorithms were compared to derive constituent absorption coefficients from 

12 remote sensing reflectance ( R rs ) ratios. The validation match-ups showed that the mean 

13 absolute percent differences (MAPD) were typically <35%, although higher errors were 

14 found for ad retrievals. Seasonal and spatial variability of satellite-derived absorption 

15 coefficients and [Chi a] was apparent and consistent with field data. CDOM is a major 

16 contributor to the bio-optical properties of the SMAB, accounting for 35-70% of total 

17 light absorption by particles plus CDOM at 443 nm, as compared to 30-45% for 

18 phytoplankton and 0-20% for non-pigmented particles. The overestimation of [Chl a] 

19 from the operational satellite algorithms may be attributed to the strong CDOM 

20 absorption in this region. River discharge is important in controlling the bio-optical 

21 environment, but cannot explain all of the regional and seasonal variability of 

22 biogeochemical constituents in the SMAB. 
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1 1 . Introduction 

2 Satellites such as the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) and the 

3 MODerate-resolution Imaging Spectrometer (MODIS) have been widely applied to the 

4 study of biogeochemical processes [IOCCG, 1999; McClain et al., 2004], Based on bio- 

5 optical theory, the satellite measurement, here remote sensing reflectance ( R rs ), is often 

6 related to inherent optical properties (IOPs) such as the absorption coefficient (a) and 

7 scattering coefficient ( b ) [Garver and Siegel, 1997; Gordon et al., 1988; Kirk, 1994; 

8 Maritorena et al, 2002; Mobley, 1994]. IOPs are often related to relevant 

9 biogeochemical constituents such as chlorophyll a concentration ([Chi a]), dissolved 

10 organic carbon (DOC), particulate organic carbon (POC), and suspended sediment 

1 1 [. Ferrari et al., 2003; Mannino et al., 2008; Rochelle-Newall and Fisher, 2002; Siegel et 

12 al., 2002], 

13 Although chlorophyll a pigment plays a critical role in understanding the bio- 

14 optical properties in oceanic waters, it is insufficient to fully characterize the 

15 biogeochemical properties, especially in coastal waters where colored dissolved organic 

16 matter (CDOM) and sedimentary resuspended matter often overwhelm phytoplankton in 

17 the contribution to bio-optical properties [Gordon and Morel, 1983; IOCCG, 2006; Kirk, 

18 1994; Mobley, 1994], In general, IOPs are composed of four components: pure water, 

19 phytoplankton, CDOM, and non-pigment ed particles [Kirk, 1994; Mobley, 1994]. 

20 Absorption from components other than pure water is often considered to be strongly 

21 correlated to [Chl a] in oceanic Case 1 waters, while such an assumption often breaks 

22 down in Case 2 waters (e.g. coastal and inland waters) [Gordon and Morel, 1983; 

23 IOCCG, 2000; Kirk, 1994; Mobley, 1994], 
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1 In theory semi-analytic (SA) models, which apply spectral deconvolution, are 

2 applicable to retrieve constituent IOPs from R rs [IOCCG, 2006], For instance, the GSM01 

3 model [ Garver and Siegel, 1997; Maritorena et al, 2002] produces [Chla], absorption 

4 coefficient of CDOM and non-pigmented particles (a^), and particulate backscattering 

5 coefficient ( bb p ). Unfortunately there are at least two problems with SA models applied to 

6 coastal waters. First, SA models require detailed knowledge of IOP relationships, which 

7 vary regionally or seasonally in coastal waters and in fact are regionally specific or 

8 empirically derived [Babin et al, 2003a, 2003b; Magnuson et al, 2004]. Second, SA 

9 models are equally sensitive to signals at all wavelengths and require them to be equally 

10 accurate. In coastal waters the satellite-derived water-leaving radiances (Z„.) at shorter 

1 1 wavelengths (e.g. 412 and 443 nm) often contain some error. Incorrect atmospheric 

12 correction due to inadequate information on aerosol absorption and the selection of 

13 inappropriate aerosol model, along with the weak signal-to-noise ratio due to strong 

14 CDOM absorption, often causes errors on L w derivation in coastal waters [Bailey and 

15 Werdell, 2006; IOCCG, 2000; Siegel et al, 2000, 2005], Consequently, the application of 

16 SA models in coastal waters faces a significant challenge because of the requirement of 

1 7 highly accurate L w at all wavelengths. 

18 Empirical algorithms do not require a full understanding of fundamental bio- 

19 optical theory. For example, they provide a direct link between satellite-sensed radiance 

20 and relevant bio-optical parameters such as [Chl a] and diffuse attenuation coefficient 

21 {Kfi) on global and regional scales [Harding et al, 2005; Mueller, 2000; O ’Reilly et al, 

22 1 998, 2000; Signorini et al , 2005], The creation of empirical algorithms, however, 

23 requires a sufficient size of highly accurate field measurements spanning all seasons and 


4 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


adequate spatial coverage for the regions of interest. Thus, empirical algorithms are 
subject to updates as the dataset increases in size. In the work presented here, a set of 
self-consistent field R n data is applied to derive absorption coefficients of oceanic 
constituents in the U.S. southern Middle Atlantic Bight (SMAB). Absorption coefficients 
are very important bio-optical properties in the study of radiative transfer modeling and 
heat budget [Mobley, 1994], carbon flux (e.g. [Chl a], primary production, DOC, and 
POC) [Arrigo and Brown, 1996; Behrenfeld et ah, 2005; IOCCG, 2006; Mannino et al, 
2008; Mavra et al, 2007; Rochelle-Newall and Fisher, 2002], water quality (e.g. diffuse 
attenuation coefficient) [Mueller, 2000], and oceanic physical processes (e.g. salinity 
distribution) [Rochelle-Newall and Fisher, 2002]. The objectives of this work were: 1) to 
develop and validate satellite algorithms in deriving constituent absorption coefficients 
for phytoplankton, non-pigmented particles, and CDOM and [Chi a] near the ocean 
surface within the SMAB to within ±35% uncertainty, 2) to determine the relative 
importance of phytoplankton, CDOM, and detritus in sunlight absorption, and 3) to 
evaluate the seasonal and regional impacts of river discharge on biogeochemical 
constituents in the SMAB. 

2. Methods 

2.1. Study region and field experiments 

This study focuses on the SMAB from the Delaware Bay (DB) mouth to the 
region south of the Chesapeake Bay (CB) mouth (Figure 1). This region is well 
recognized for the significant impacts by riverine discharge from the Delaware and 
Chesapeake Bays, which account for most of the salinity variability of the SMAB [Acker 
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1 et al. , 2005; Austin, 2002; Harding, 1994]. The magnitude of freshwater run-off, along 

2 with wind and tidal forcing, generates periodic outflow plumes (e.g. winter-spring plume 

3 and fall sub-plume) for this region, and enhances the bio-optical complexity of the 

4 SMAB compared to pelagic regions of the Atlantic Ocean [Acker et al., 2005; Harding, 

5 1 994; Johnson, 200 1 ; O ’Reilly and Zetlin, 1998; Rennie et al., 1999] . 

6 Multiple cruises were conducted in this region during 2004-2006, including the 

7 Bio-physical Interactions in Ocean Margin Ecosystems cruises (BIOME) during 30 

8 March to 1 April 2005 (BIOME1), 26 to 30 July 2005 (BIOME2), 9 to 12 May 2006 

9 (BIOME3), and 2 to 6 July 2006 (BIOME4), and the Chesapeake Bay Plume cruises 

10 (CBP) during 27 May and 3 November 2005, and 6 September and 28 November 2006, 

1 1 and the Chesapeake Bay Hydrological survey (CBH) during 5 May, 5 July, 1 September, 

12 15 October, and 15 November 2004, and 10 January, 26 May, 21 June, 19 August, and 23 

13 September 2005 (Figure 1). The collected bio-optical data included, but not limited to, 

14 phytoplankton pigments, IOPs (absorption), and R rs (only on BIOME and CBP cruises in 

15 2005). Water samples were collected at multiple depths with Niskin bottles. 

16 

17 2.2. Pigments and absorption coefficient 

18 Pigment samples were collected on 25 mm GF/F filters under a gentle vacuum 

19 (<5 in Hg) and stored in liquid nitrogen in the field before transfer to a -80°C freezer in 

20 the laboratory. Pigments were analyzed at Horn Point Laboratory by reverse-phase high- 

21 performance liquid chromatography (HPLC) with a C8 column on the HPLC system 

22 equipped with photodiode array detector [Van Heukelem and Thomas, 2001], [Chl a] 
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1 was calculated as the sum of concentration from monovinyl Chla, divinyl Chla, and 

2 chlorophyllide a. 

3 Absorption coefficients of particles ( a p ) and non-pigmented particles (ad) were 

4 determined by the quantitative “filter pad” method following the recommendations of 

5 Mitchell et al. [2002], Particulate samples were collected on 25 mm GF/F filters under a 

6 gentle vacuum (<5 in Hg) and stored in liquid nitrogen in the field before transfer to a - 

7 80°C freezer in the laboratory. Non-pigmented particulate samples were defined as the 

8 detritus component of particulate samples after two cold methanol extractions (first 5 ml 

9 for 10 minutes, then 10 ml for 1 hour) [Kishino et al, 1985], Artificial seawater pre- 

10 filtered through 0.2 pm Whatman Nuclepore filters was used to rinse off methanol and to 

1 1 hydrate the GF/F filters of both blanks and samples. CDOM samples were collected by 

12 filtering seawater through pre-combusted (6 hours at 450°C) GF/F filters and stored 

13 under refrigeration (4 to 8 °C). In the laboratory, CDOM samples were warmed to room 

14 temperature and filtered through 0.2 pm Whatman Nuclepore (polycarbonate) or Gelman 

15 Supor (polyethersulfone) filters prior to analysis [Mannino et al, 2008], 

16 Absorbance spectra were measured using a double-beam Cary 100 Bio 

17 Ultraviolet -Visible scanning spectrophotometer through 250-800 nm (CDOM in Suprasil 

18 quartz 10 cm pathlength cells) or 300-800 nm (particles) in 1 nm intervals. Blank GF/F 

19 filters hydrated with 0.2 pm pre-filtered artificial seawater and ultraviolet (UV) oxidized 

20 Milli-Q water were used as the blank and reference for particulate absorbance and 

21 CDOM absorbance, respectively. Null correction was made by subtracting the mean of 

22 absorbance at 790 to 800 nm for particulate samples for each spectrum. No null 

23 correction for CDOM absorption was made since the raw absorbance of seawater samples 
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1 at 690-700 nm was within the noise level of the instrument [. Mannino et al., 2008], The 

2 multiple-scattering effect for calculating particulate absorption coefficient ( a p ) and non- 

3 pigmented particulate absorption coefficient {af) was corrected following the method of 

4 Mitchell [1990], from which non-pigmented particulate samples were assumed to have 

5 the same multiple-scattering amplification factor to total particulate samples [Mitchell et 

6 al. , 2002], Phytoplankton absorption coefficient (a p i,) was calculated as a ph = a p -a d . 

7 The absorption coefficient by CDOM and non-pigmented particles ( ad g ) was calculated as 

8 the sum of ad and CDOM absorption coefficient ( a g ). Total absorption coefficient (a) was 

9 calculated as a = a w +a p +a g , where pure water absorption coefficient (n„.) was adopted 

1 0 from Pope and Fry [1997], 

1 1 The absorption coefficient from non-pigmented particles, CDOM, or their sum 

12 ( a x ) was fitted to an exponential function as: 

13 

14 a x (A) = a x (A 0 )exp[-S x (A-Jl 0 )] (1) 

15 

16 Here, S x represents the exponential slope for absorption coefficient from non-pigmented 

17 particles ( Sd ), CDOM (S g ), or their sum (Sd g ). We selected the reference wavelength, ko, 

18 equal to 380 nm. In Equation (1), the wavelengths analyzed were 350 to 600 nm for a g {X), 

19 and 380 to 730 nm for ad{k) but excluding 400 to 480 nm and 620 to 710 nm to avoid the 

20 chlorophyll pigment peaks due to methanol’s incapability to extract some pigments as 

21 discussed by Jeffrey et al. [1997], and 380 to 600 nm (but excluding 400 to 480 nm) for 

22 ad g {k ) [Babin et al., 2003b]. 
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1 The phytoplankton absorption coefficient ( a p t , ) is typically related to [Chl a] as a 

2 power function [Bricaud etal., 1995, 1998; Prieur and Sathyendranath, 1981]: 

3 

4 a ph (A) = A 0 (A)[Chl_a] A ' w (2) 

5 

6 The modification of Equation (2) provides an expression of a p h from its value at a 

7 reference wavelength (here 670 mn): 

8 

9 a ph (A) = B 0 (A)[a ph (670)] B ' w (3) 

10 

1 1 Ai(A.) and Bj(A.) are derived coefficients. [Chi a] can also be determined from a p i,(670) by 

12 a power function similar to Equation (3). 

13 

14 2.3. Apparent optical properties from in situ measurements 

15 The remote sensing reflectance ( R rs ) spectra (bands centered at 320, 340, 380, 

16 395, 412, 443, 465, 490, 510, 532, 555, 560, 625, 665, 670, 683, 710, 780, and 860 nm, 

17 and each band is 10 nm wide at Hill-width half max) were determined with a BioPro in- 

18 water profiling spectroradiometer (Biospherical Instruments, Inc., San Diego, CA), as 

19 described in detail by Mannino et al. [2008], The instrument was deployed multiple times 

20 for each station, and the absolute uncertainty was less than 5%. The R rs at 55 1 mn was 

21 calculated from a linear interpolation of values at 532, 555, and 560 nm [Mannino et al., 

22 2008], The R rs at 488 nm was assumed to be equivalent to the value at 490 nm. R rs 

23 measured at 6 stations during BIOME 1 cruise, 1 9 stations during BIOME2 cruise, and 3 


9 



1 and 6 stations during CBP cruises on 27 May and 3 November 2005 were included for 

2 analysis in this paper. 

3 

4 2.4. Satellite ocean color validation 

5 The method to process satellite images was described in Mannino et al. [2008] 

6 following Bailey and Werdell [2006] protocols. SeaWiFS and MODIS-Aqua observations 

7 were processed from Level 1 to Level 2 using the SeaWiFS Data Analysis System 

8 software (SeaDAS version 5.1.1 and msll2 version 5.6.3). The pixels were masked after 

9 atmospheric correction by any of the following flags: land, cloud or ice, high top-of- 

10 atmosphere radiance, low normalized water-leaving radiance at 55 1 or 555 nm, stray 

1 1 light, sun glint, or atmospheric correction failure) [. Bailey and Werdell, 2006]. Pixels with 

12 L„.(4 12)<0.2 mW cm' pm’ sr' were excluded to minimize the impacts from atmospheric 

13 over-correction in causing negative or significantly reduced water-leaving radiance 

14 [ Siegel et al, 2002]. The 3x3pixel arrays centered on the field stations, each with ~1 km 

15 resolution (sensor native), were analyzed. The satellite observations which occurred 

16 within ±8 hours and ±32 hours of in situ measurements were considered for match-up 

1 7 analysis to allow for inclusion of sufficient data points. 

18 

19 2.5. Analysis and validation methods 

20 Three curve-fitting functions, linear function (Model II regression) after log- 

21 transformation (loglinearmodel; Equation 4), fourth-order polynomial function after 

22 log-transformation similar to OC4V4 \0 ’Reilly et al, 1998, 2000] 

23 (log_polynomial_model; Equation 5), and one-phase exponential decay function 
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1 (exponentialmodel; Equation 6), were developed to correlate R rs band ratio to the 

2 relevant absorption products: 

3 

4 log[a,.(/l)] = C 0 (/l) + C 1 (/l)^ (4) 

5 

6 log[a,(/i)] = D 0 (A) + D X (A)R + D 2 (A)R 2 + D 3 (A)R 3 + D 4 (A)R 4 (5) 

7 

8 a,(X) = G 0 (A) + G 1 (A)exp[-G 2 (A)4=^4] (6) 

Z3 (4-2 ) 

9 

10 Here, R = log[i? re (A x ) / R n (A 2 )] , and Ay and A? represent the various bands evaluated, and 

1 1 C„ Dj and G, are wavelength-specific derived coefficients, and a, is the analyzed 

12 absorption coefficient as a p h, o<y, a g , or a<\ g . In addition, a one-phase exponential function 

13 to determine R rs from a, similar to Mannino et al. [2008] was also developed, and a, was 

14 then calculated by its reverse function (reverse exponential model; Equation 7): 

15 

1 6 4=777 = HM) + H x (A) exp [-ff 2 (A)a, (A)] (7) 

% (4-2 ) 

17 

18 Four products, ^,(670), 0^380), <7^(380), and <7^380), whose surface measurements 

19 were represented as the site values, were analyzed with the above equations. The mean 

20 absolute percent difference (MAPD) and root mean square error (RMSE) between the 

21 modeled products (C a i g ) and field measurements ( Cf ie id ) were calculated. 
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2 MAPD = ^ (Cn ' g x 1 00% (8) 

N 

3 

4 *MKE = V2X. t-Cj^f/N (9) 

5 

6 2.6. Monthly time series analysis 

7 Monthly Level 3 mapped MODIS-Aqua images (4 km resolution) from July 2002 

8 to December 2006 were downloaded from the NASA ocean color website 

9 (http://oceancolor.gsfc.nasa.gov ) on 8 January 2008. Empirical algorithms developed 

10 from this paper were applied to calculate products such as [Chla], a p h, a^, and a g . Three 

1 1 stations— Location A (75.90W, 36.93N), B (75.30W, 36.93N), and C (74.77W, 36.93N) 

12 representing a transect from the Chesapeake Bay mouth to an outer shelf location — were 

13 selected and plotted to demonstrate a monthly time series. 

14 

15 3. Results and Discussion 

16 3.1. Absorption spectra 

17 We observed seasonal transitions in phytoplankton absorption coefficients in 

18 April-May and October-November periods. Therefore, at least two seasonal algorithms 

19 (May-October and November- April) are required to describe phytoplankton absorption 

20 relationships (Figure 2a and Table 1). In general, the phytoplankton absorption ratio 

21 [a p h(P)/a p i, (670)] in May-October is higher than that in November- April (Figure 2a). This 

22 coincides with historical observations that the dominant phytoplankton taxa are diatoms 
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1 in winter and spring, but the phytoplankton assemblage transitions to a greater proportion 

2 of dinoflagellates, cryptophytes, and cyanobacteria in summer and fall [Adolf et al, 2006; 

3 Marshall and Alden, 1993], Although phytoplankton absorption coefficients are subject 

4 to seasonal variation, the relationship between a p h(670) and [Chla] (r =0.964, N=230; 

5 Figure 2b) is relatively constant seasonally due to the dominant contribution from Chl a 

6 to a p h(6 70) [Jeffrey et al, 1997]: 

7 

8 [Chi _ a] = 70.632 x [a ph (670)] 1 184 (10) 

9 

10 Seasonal variability of riverine discharge, along with other factors such as wind forcing 

1 1 and direction, may cause the seasonal variability of phytoplankton taxonomic 

12 composition, pigment package effect, and therefore normalized phytoplankton absorption 

13 spectra (Figure 2a) [Babin et al, 2003b; Bricaud et al, 1995, 1998; Trees et al, 2000], 

14 The pigment package effect refers to a consequence of the fact that in the natural waters 

1 5 pigment molecules are not uniformly distributed but are contained within discrete 

16 packages such as chlorop lasts, cells, and cell colonies, which causes a flattening of the 

17 phytoplankton absorption peak due to self-shading wherever pigments are localized 

18 within cell membranes [Duysens, 1956]. 

19 Equation (1) described the exponential decay characteristics of a d , a g , and a dg 

20 rather well with coefficients of determination (r ) of >0.95 for ad, and >0.99 for a g and a dg . 

21 The exponential slopes ( S) covered a wide range with mean ± 1 standard deviation of 

22 0.0122±0.0023 (N=247, ranging from 0.0084 to 0.0260) for a d , 0.0170±0.001 1 (N=300, 

23 ranging from 0.0137 to 0.0221) for a g , and 0.0148±0.0014 (N=222, ranging from 0.0122 
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1 to 0.0205) for ad g . However, the general exponential decay relationships of S (ordinate) 

2 versus absorption coefficients (abscissa) for coastal regions that are significantly 

3 impacted by freshwater discharge as suggested by Carder et al. [1989] was not apparent 

4 for the SMAB dataset. Part of the reason was that the dataset presented in this paper 

5 (Figure 1) did not extend far into the estuaries to cover the full range of water types from 

6 freshwater to oceanic waters. The observations above suggest that using mean S values 

7 may cause significant errors in predicting a x spectra. Multiple algorithms to retrieve ad, a g , 

8 and ad g at multiple wavelengths may be required, as discussed in the next section. 

9 

10 3.2. Algorithm development and validation 

11 Three R rs band ratios (412/555, 443/555, and 490/555 for SeaWiFS, and 412/551, 

12 443/551, and 488/551 for MODIS-Aqua) were compared to determine the best algorithm 

13 performance. No matter which of the curve-fitting models from Equations (4)-(7) was 

14 selected, algorithms based on R rs band ratio at 490/555 or 488/55 1 performed similar to, 

15 or better than, the other two band ratios. Figure 3 shows examples of model performance 

16 on predicting a p h(610), 0^(380), ag(380), and aj^, 380) from Equations (4)-(7), 

17 respectively. Performance from other models, i.e. Equation (4) on predicting a</(380), 

18 a g (380), and av g (380), showed similar results to Figure 3 (data not shown). Due to the 

19 possibly poor quality of satellite water-leaving radiance at shorter wavelengths (e.g. 412 

20 and 443 nm) in coastal waters [ Bailey and Werdell, 2006; Siegel et al, 2002], we selected 

21 R rs band ratios at 490/555 or 488/551 for further analysis. Although the selection of only 

22 two bands will cause interdependence of satellite-derived products, it still provides 
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1 valuable information on bio-optical properties of the SMAB in the absence of appropriate 

2 semi-analytic algorithms based on more bands. 

3 The loglinearmodel (Equation 4) proved suitable to validate satellite-derived 

4 a p h, ad, and cid g with relatively high r and relatively low MAPD and RMSE (Figures 4-5 

5 and Table 2). The r 2 for the log linear model ranged from 0.92 to 0.93 for all three 

6 products of a p i,(670), arf(380), and 0^(380), while 0.84 to 0.95 for the other three models. 

7 The MAPD for the log linear model was typically similar to, or lower than, the other 

8 three models, regardless of the satellite sensor (SeaWiFS or MODIS-Aqua) and overpass 

9 satellite/in situ time window selected (±8 hours or ±32 hours) (Figures 4a-b and 5a-b). 

10 This was also supported by the validation results from the RMSE comparisons (Figures 

1 1 4c-d and 5c-d). The exponential model may provide lower MAPD and RMSE in some 

12 cases (Figures 4-5), but was not selected due to its relatively low r (0.84 to 0.87). In 

13 contrast, the reverse exponential model yielded the best validation results for ag(380) 

14 with the highest r 2 (0.90 versus 0.77-0.86) and typically lowest MAPD and RMSE 

15 (Figures 4-5 and Table 4) as compared to the other models. Mannino et al. [2008] 

16 showed the same model for a g but did not include the CBH stations in their validation 

17 analysis. The selected R rs band ratio models can also be applied to derive ad, a g , and ad g at 

18 multiple wavelengths (Tables 2-3). Except for the higher MAPD for ad (34.8-57.5% for 

19 SeaWiFS and 41 .9-65.3% for MODIS-Aqua), the selected regression methods typically 

20 limited MAPD for a g and ad g to within 30% for wavelengths between 350 and 555 nm 

21 (Table 4). The exponential decay slopes ( S) for ad, a g , and ad g can be derived from non- 

22 linear regression methods with R rs band ratio models at multiple wavelengths (e.g. 355, 

23 380, 400, 412, 443, 490, 510, 531, and 555), and agree reasonably well with field 
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1 derivations (Table 4). The GSM01 [Garver and Siegel , 1997; Maritorena el al., 2002] 

2 and its regional version (GSM01-CB) [Magnuson et al, 2004] resulted in relatively high 

3 MAPD and RMSE when compared with the field measurements (Figures 4-5). The 

4 GSM01 model was developed for global ocean application, and thus is not optimized for 

5 the variability of in-water constituents observed in near-shore coastal regions, such as 

6 variable or region-specific S values. However, even the optimized GSM01 model for this 

7 coastal region (GSM01-CB; primarily Chesapeake Bay and near-shore coastal ocean) 

8 developed by Magnuson et al. [2004] did not perform significantly better than GSM0 1 , 

9 e.g. MAPD=42.6% and 46.0% for o^(670) and a^(380) from GSM01-CB, versus 43.3% 

10 and 40.7% from GSM01 for SeaWiFS ±8 hour overpass window (Figures 4-5). The lack 

11 of adequate knowledge of backscattering coefficients, as well as the higher uncertainty of 

12 R rs at shorter wavelengths (e.g. 412 and 443 nm) from satellite measurements, may 

13 account for the performance of GSM01 and GSM01-CB for this region [Bailey and 

14 Werdell, 2006; IOCCG, 2000; Magnuson et al, 2004; Siegel et al, 2000, 2005]. 

15 By applying the regression results shown in Table 1, phytoplankton absorption 

16 coefficients at other visible wavelengths can also be derived from a p i,(6 70), which can be 

17 derived from satellite radiance observations as shown in Table 2. The validation match- 

18 ups based on this approach yielded similar MAPD accuracy levels for a p h{X) at 412, 443, 

19 488, 490, 510, 667, and 678 nm to o^/,(670) (21.5-26.1% versus 25.7% for SeaWiFS and 

20 19.0-28.1% versus 21.2% for MODIS-Aqua), but relatively higher MAPD at 531, 551, 

21 and 555 nm (27.2-30.9% for SeaWiFS and 33.1-43.1% for MODIS-Aqua) due to the 

22 relatively higher measurement errors from weaker absorption at these wavelengths (Table 

23 4). If we exclude those stations with extremely low a p h (e.g. <0.003 m' 1 at 555 nm), the 
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1 MAPD for a p h(k ) was at the same accuracy level for all visible wavelengths evaluated. 

2 The satellite derivation of [Chla] from Equation (10) had slightly lower but still 

3 reasonable accuracy level compared to that for a p i, (67 0) with MAPD of 32.3±28.2% 

4 (N=29) for SeaWiFS and 28.8±20.6% (N=14) for MODIS-Aqua (Table 4). 

5 The validation match-ups between field measurements of absorption constituents 

6 and satellite derivations within ±8 hours demonstrated reasonable agreement (Figures 6- 

7 7). The SeaWiFS and MODIS-Aqua match-ups have similar r 2 varying from 0.80-0.97, 

8 and show slope values from 0.52 to 0.73 and from 0.46 to 0.80, respectively. Including 

9 the stations applied to develop the algorithms, which increases dataset size by 10-12 

10 points for SeaWiFS and 5-6 points for MODIS-Aqua, yielded improvements in all the 

11 slopes (e.g. 0.73-0.87 for SeaWiFS and 0.57-0.64 for MODIS-Aqua). When extending 

12 the match-up dataset from ±8 hours to ±32 hours of the satellite overpass window, similar 

13 improvements were also found due to significant increase of dataset size by 50% to 140% 

14 (data not shown). The validation match-ups between field measurements of [Chl a] and 

15 satellite derivations after log-transformation agreed well with r =0.86 to 0.95, slope=0.87 

16 to 0.96, and RMSE=0.20 to 0.24 (Figure 7). Similar statistical results were found for 

17 a p h( 670) when including those stations used to develop the algorithms (data not shown). 

18 It implies that these empirical algorithms are relatively successful and should improve as 

19 the size of the dataset increases. 

20 Above all, the satellite-derived absorption coefficients from selected functions 

21 (loglinearmodel for a p h, and ajg, and reverseexponentialmodel for a g ) yielded 

22 relatively good results for the SMAB. Since the dataset used for algorithm development 

23 did not include stations from the CBH cruises (because AOP data was not collected at 
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1 those stations), the validation analyses should improve after excluding those CBH 

2 stations. For example, MAPD improved from 25.7% to 21 .9% for a p h(6 70), from 35.3% 

3 to 29.5% for ad( 380), from 25.2% to 20.9% for Og(380), from 26.5% to 23.8% for 

4 dd g { 380), and from 32.3% to 24.2% for [Chl a] for SeaWiFS ±8 hours overpass time 

5 window. Therefore, in the future the addition of complete datasets from the lower bay 

6 locations will expand the dynamic range of the algorithms and may significantly improve 

7 the model capability in more turbid areas of the SMAB. 

8 

9 3.3. Seasonal variability 

10 These empirical algorithms can be applied to study the spatial and seasonal 

1 1 variability of coastal ocean constituents. Figure 8 shows typical examples of the spatial 

12 distribution of the absorption coefficients within the SMAB during four seasons: summer 

13 (June-August), fall (September-November), winter (December-February), and spring 

14 (March-May). The satellite images clearly show the gradients from high to low 

1 5 constituent concentration between the coast and the open ocean as well as the 

16 riverine/esturine outflow impact along the coast (Figure 8). The seasonal variability of 

17 phytoplankton absorption [e.g. a p h{ 443)] may be due primarily to the river discharge rate 

18 from the bay mouths [Acker et al., 2005; Adolf et al., 2006; Marshall and Alden, 1993; 

19 Marshall et al. , 2006], The monthly mean flow rates out of the Chesapeake Bay for these 

20 selected images were 510, 1648, 2983, and 1463 m 3 s' 1 for August 2005, November 2005, 

21 February 2006, and May 2006, respectively (Data sources: 

22 http://waterdata.usgs.gov/nwis/; written communication from Gary Fisher, U.S. 

23 Geological Survey, 17 July 2007). Consequently, lower phytoplankton abundance 


18 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


occurred during the dry season in summer 2005 compared to the other three seasons. The 
seasonal variability of detritus absorption [e.g. 0^443)] is complex. At least two primary 
sources of detritus from riverine/esturine outflow and sedimentary resuspension 
controlled ad in the SMAB. The significant contribution from storm-driven sedimentary 
resuspension in autumn may explain the higher detritus content in November 2005. The 
seasonal variability of CDOM absorption [e.g. Og( 443)] may be controlled primarily by 
the degree of riverine inputs of degraded terrestrial vegetation to the SMAB [Del Vecchio 
and B lough, 2004; Mannino et al., 2008]. 

CDOM plays a critical role in contributing to sunlight absorption and thus impacts 
primary production in the SMAB by reducing the amount of photosynthetically active 
radiation (PAR) available for phytoplankton growth [. Air i go and Brown , 1996], At 443 
nm, CDOM accounted for 35-70% of total light absorption excluding water absorption 
(a pg ), as compared to 0-20% for non-pigmented particles, and 30-45% for phytoplankton 
(Figures 9-10). Pure seawater absorption [a„.(443)~0.007 m' 1 ; Pope and Fry, 1997] 
typically accounts for a negligible fraction (~3.0%) of o(443) in the SMAB. The 
relatively low contribution of detritus absorption within coastal ocean regions was also 
reported by Siegel et al. [2002] and may explain the performance of the satellite-derived 
ad(k) in the validation analysis (Table 4, and Figures 4-6). The satellite-derived adQC) 
underestimates ad(k) with respect to field samples collected in the near-shore ocean 
region (<20 m bottom depth) and overestimates a ( i(k) in water with very low aj(k) 
(Figures 6, 8, and 10). The gradients of high to low percentages of a p h and ad, and low to 
high percentage of a g from the coast to the open ocean were consistent with field 
measurements (Figure 10). During the dry season (e.g. August 2005), CDOM accounts 
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1 for a higher percentage of total absorption than during the wet season (e.g. February 2006) 

2 (Figure 9). Such a phenomenon may be explained by the impact of river discharge in 

3 contributing nutrients to support phytoplankton growth as well as the export of terrestrial 

4 CDOM. During the wet season, CDOM and phytoplankton abundance are both elevated, 

5 but phytoplankton blooms increase the relative percentage of phytoplankton absorption 

6 compared to CDOM absorption. During the dry season, CDOM and phytoplankton are 

7 both low, but the reported higher primary productivity and mature grazer community may 

8 result in a higher percentage of phytoplankton to be grazed and degraded which in turn 

9 reduces the relative percentage of phytoplankton absorption [Adolf et al, 2006; Marshall 

10 and Nesius, 1996; Marshall et al, 2006] and increase CDOM through grazer and 

1 1 microbial processing of organic matter [Nelson et al, 2004; Steinberg et al, 2004], Since 

12 phytoplankton pigments have a much weaker relationship with CDOM absorption than 

13 with phytoplankton absorption, the significant contribution of CDOM absorption may 

14 pose complications for applying global operational algorithms (e.g. OC4V4 and OC3M) 

15 [O 'Reilly et al, 1998, 2000] to coastal regions. For CDOM -rich Case 2 waters such as 

16 the Chesapeake Bay, OC4V4 has been found to significantly overestimate [Chla], 

17 especially for offshore regions of the SMAB [Harding et al, 2005; Magnuson et al, 

18 2004], Our results also support this conclusion. For example, match-ups within ±8 hours 

19 showed that OC4V4 performed better for the lower CB region (e.g. CBH stations) with 

20 MAPD of 33.8% as compared to 79.6% for whole SMAB region (data not shown). It 

2 1 implies that the relative difference between our approach and operational algorithms 

22 would be relatively small in near-shore regions but high in offshore regions. The spatial 

23 distribution and the seasonal variability of [Chl a] based on our approach displayed 
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1 similar trends as those from OC4V4 and OC3M algorithms, but significantly reduced the 

2 overestimation by operational Chla algorithms in the offshore region of the SMAB 

3 (Figure 1 1). In general, the ratios of [Chl a] based on operational Chl a algorithms to our 

4 approach increase with the increase of CDOM contribution to light absorption (Figures 9 

5 and 11). OC4V4 and OC3M [Chl a] were higher by 0-0.5 times for the inner-shelf 

6 region, 0.4- 1.2 times for the middle shelf region, and 1-2 times for the outer shelf region 

7 (Figure 1 1). This higher ratio trend toward offshore demonstrates the impact of CDOM 

8 on ocean color products in the SMAB. 

9 The satellite derivations of absorption coefficients provide tools to study 

10 biogeochemical processes and radiative transfer. For examples, DOC and salinity can be 

1 1 strongly correlated to CDOM absorption [Del Vecchio and Blough, 2004; Mannino et al, 

12 2008; Rochelle-Newall and Fisher, 2002], and primary productivity is correlated to 

13 phytoplankton absorption [Behrenfeld et al, 2005; Marra et al, 2007], The knowledge of 

14 absorption also provides methods to study other IOPs from space. For example, the 

15 expression of R rs from absorption and backscattering [Garver and Siegel, 1997; Gordon 

16 et al, 1988; Maritorena et al, 2002] and the empirical expression of absorption from R rs 

17 ratio make it possible to express backscattering into R rs . The knowledge of backscattering 

18 might significantly improve the capability of semi-analytical models in deriving ocean 

19 color products from space [Magnuson et al, 2004]. 

20 Although we have shown the significant impact of river discharge on 

21 biogeochemical constituents in the SMAB, the direct link between them should be 

22 interpreted with caution. First, the impact of river discharge on the coastal region of 

23 Chesapeake Bay is different from Delaware Bay. The lower Chesapeake Bay is subject to 
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1 nutrient limitation for phytoplankton growth, in contrast to light availability in the 

2 Delaware Bay [Harding et al, 1986; Marshall and Alden, 1993], Therefore, an increase 

3 in river discharge is more likely to cause a phytoplankton bloom in the lower Chesapeake 

4 Bay by driving more nutrients downstream, while an increase in turbidity from higher 

5 river discharge may decrease primary production in the lower Delaware Bay. Second, the 

6 impact of river discharge is subject to seasonal variability and distance from the bay 

7 mouths, as shown in the following for the coastal region of Chesapeake Bay. In the inner- 

8 shelf region, the correlation coefficients (r) between river discharge rate and 

9 biogeochemical products are low (e.g. r=0.05-0. 12 for [Chi a] and -0.01-0. 14 for a g ) for 

10 all seasons except for summer (Figure 12). The poor correlation may be due to averaging 

1 1 out the higher frequency responses (less than one week) for the export of nutrients and 

12 CDOM, respectively, from the bays. During summer, the vertical stratification is well 

13 developed [ Verity et al, 2002], and the strength of river discharge represents the flux of 

14 nutrients for phytoplankton growth. Thus, the correlation between biogeochemical 

15 products and river discharge improves for summer (e.g. r=0.48 for [Chi a] and a g ) 

16 (Figure 12). In the middle shelf region, however, river discharge is significantly 

1 7 correlated to biogeochemical products in winter but poorly correlated during other 

18 seasons (e.g. r=0.79 for [Chl a] and a g in winter and r=0. 09-0.36 in other seasons) 

19 (Figure 12). During winter, low water temperature and a less mature grazer community 

20 may cause the phytoplankton biomass to be linked directly to nutrient availability, which 

21 is driven primarily by river discharge and by wind-induced vertical mixing of nutrients 

22 from depth [Adolf et al., 2006; Marshall and Alden, 1993]. As the zooplankton and 

23 bacterial communities develop into spring and summer, lower phytoplankton biomass 
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1 and higher primary productivity are expected [Adolf et a!., 2006] and the direct response 

2 of the biological system to river discharge dissipates. The outer shelf region shows a 

3 similar pattern but a lower correlation coefficient (e.g. r=0.49 for [Chl a] and a g in 

4 winter, and -0.13-0.31 in other seasons) with river discharge than the middle shelf region 

5 (data not shown). 

6 Other physical factors than river discharge, such as water temperature and wind 

7 forcing, anthropogenic activities, and even climate change, can also impact 

8 phytoplankton abundance, productivity, and carbon distributions in the SMAB. For 

9 example, the direction and distribution of the Chesapeake Bay plume is highly dependent 

10 on the wind stress direction. During winter and early spring northerly winds 

1 1 (downwelling favorable) and the along-shore southward current force the Chesapeake 

12 Bay [Rennie et ai, 1999; Verity et ai, 2002] and Delaware Bay [Sanders and Gan’ine, 

13 2001] plumes to flow southward along the coast. As winds reverse later in spring the 

14 southerly along-shore flow weakens, and the Chesapeake Bay plume broadens and flows 

15 offshore, primarily to the south and east. Upwelling-favorable conditions can initiate 

16 local phytoplankton blooms and contribute additional particles to surface waters [Johnson 

17 et al. , 2001]. The Chesapeake Bay estuarine ecosystem has experienced a large increase 

18 in anthropogenic nutrient loading and reductions in the past half century which have 

19 affected the floral composition and biomass [Harding, 1994; Paerl et ai, 2006]. 

20 Furthermore, climate forcing (e.g. hurricanes, drought, etc.) significantly influences 

21 phytoplankton dynamics (e.g. by reducing vertical stratification, increasing sedimentary 

22 resuspension, and redistributing particles from hurricane forcing) [Miller and Harding, 

23 2005; Paerl et ai, 2006]. 
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1 Above all, the impacts from physical factors (e.g. river discharge, wind forcing, 

2 and bathymetry) on bio-optical constituents (e.g. [Chla] and a g ) are complicated and 

3 cannot be explained by a single factor [Harding, 1994], Nevertheless, we found that the 

4 variability of an optical property, the diffuse attenuation coefficient at 490 nm (K490), 

5 represents the variability of multiple bio-optical constituents. In the Chesapeake Bay 

6 inner-shelf site, the correlation coefficient (r) of [Chi a], a p h, ad, a g , and aj g to K490 was 

7 0.66-0.74, while 0.94-0.99 in the middle shelf location, and 0.93-0.98 in the outer shelf 

8 locations. These results also imply that absorption is the dominant contributor to the 

9 diffuse attenuation coefficient at offshore locations but scattering contributes 

10 significantly at near-shore locations. 

11 

12 4. Conclusions 

13 Several important conclusions can be made from the present analyses of 

14 absorption coefficients and [Chl a] derived from ocean color remote sensing. The 

15 empirical algorithms demonstrate successful retrieval of absorption coefficients and 

16 [Chl a] within a reasonable uncertainty (e.g. ±35%), and demonstrate significant 

17 improvements from the standard semi-analytic model (e.g. GSM01 and GSM01-CB) and 

18 operational algorithms (e.g. OC4V4 and OC3M). Field observations and satellite 

19 derivations both demonstrate that CDOM is the major contributor to water-column light 

20 absorption at shorter wavelengths (e.g. <500 nm), especially during the dry seasons and 

21 on the outer shelf where it can account for 35-70% of absorption by particles plus CDOM 

22 at 443 nm. River discharge plays a principal role in controlling the distribution of 

23 biogeochemical constituents, but is subject to seasonal and regional variability. 
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1 Figure captions 

2 

3 Figure 1. Map of the study area within the southern Middle Atlantic Bight (SMAB). 

4 Symbols representing the sampling stations from the following cruises are: A - BIOME1 

5 (30 March to 1 April 2005), V - BIOME2 (26 to 30 July 2005), O - BIOME3 (5 to 9 

6 May 2006), X - BIOME4 (2 to 6 July 2006), □ - CBP (four daily cruises), and ■ - 

7 CBH (ten daily cruises). 

8 

9 Figure 2. Log-transformed linear regression of phytoplankton absorption coefficient at 

10 670 nm [a / ,/,(670)] to a) a^/,(443), and b) chlorophyll a concentration [Chla]. The solid 

1 1 and dashed lines in the upper figure represent the regression for summer-fall season (May 

12 to October) and winter-spring season (November to April), respectively. 

13 

14 Figure 3. Absorption algorithms derived from field observations of remote sensing 

15 reflectance ( R rs ) from a) log linear model for a p i,( 670): log [a ph (A)] = C 0 (A) + C X (A)R , 

16 where R = log[i? rj (T, ) / R rs (A 2 )] ; and b) log_polynomial model for non-pigmented 

17 particulate absorption coefficient at 380 nm [tf,X380)]: 

18 log[a d (A)\ = D 0 (A) + D l (A)R + D 2 (A)R 2 + D 3 (A)R 3 + D 4 (A)R 4 ; and c) 

19 exponentialmodel for CDOM absorption coefficient at 380 nm [o g (380)]: 

20 a g (A) = G 0 (A) + G x {A) exp[-(j 2 (A) — — * ] ; and d) reverse exponential model for 

R n(* 2 ) 

21 absorption coefficient by non-pigmented particles plus CDOM at 380 nm [^(380)] 
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R (/L ) 

1 — - ^ 1 = H 0 (T) + H l (/l) exp [-//, (A)a dg (T)] • Regression lines from R rs band ratios of 

R ,s (4z ) 

2 412/555, 443/555, and 490/555 are represented as solid, dotted, and dashed lines, 

3 respectively. 

4 

5 Figure 4. Validation results comparing SeaWiFS observations with field measurements of 

6 a p },{ 670), a</(380), a^(380), and 380) from multiple models (loglinearmodel, 

7 log_polynomial model, exponentialmodel, reverse exponential model). Figures a) and 

8 b) show the mean absolute percent difference (MAPD), while c) and d) show the root 

9 mean square error (RMSE) of the validation results within 8 hours and 32 hours of the 

10 satellite overpass, respectively. The data from stations applied to develop the algorithms 

1 1 were not included in this analysis. The satellite derived a p h(6 70) and 0^(380) from 

12 GSM01 model [. Maritorena et al., 2002] and GSM01-CB model [Magnuson et al, 2004] 

13 are also shown for comparison. 

14 

15 Figure 5. Validation results comparing MODIS-Aqua observations with field 

16 measurements of <3^(670), <3^(380), «g(380), and <3^(380) from multiple models. See 

17 Figure 4 for details. 

18 

19 Figure 6. Comparisons of SeaWiFS and MODIS-Aqua and field observations of the 

20 absorption coefficients of a) a p h{61Qi), b) <3^(380), c ) a g (380), and d) a^g(380). The values 

21 are plotted on log scale. The satellite derivations of <3^(670), <3^(380), and <3^(380) were 

22 from the log linear model, while <3^(380) from the reverse exponential model. The 

23 match-ups procedure is limited to within ±8 hours, and the data from stations used to 
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1 develop algorithms are excluded for validation analyses. The statistical results are based 

2 on log-transformation of the data and shown on upper left for SeaWiFS and lower right 

3 for MODIS-Aqua. The solid lines represent the 1:1 lines, while dashed lines and dotted 

4 lines represent the regression for SeaWiFS and MODIS-Aqua respectively. 

5 

6 Figure 7. Comparisons of SeaWiFS and MODIS-Aqua and field observations of [Chi a] 

7 for satellite overpass window of a) ±8 hours and b) ±32 hours. The data from stations 

8 used to develop algorithms of a p h are excluded for this analysis. The solid lines represent 

9 the 1 : 1 lines, while dashed lines and dotted lines represent the regression for SeaWiFS 

10 and MODIS-Aqua respectively. See Figure 6 for detail. 

11 

12 Figure 8. The distribution of a p h, ad, and a g at 443 mn within the SMAB for 5 August and 

13 3 November 2005, and 15 February and 12 May 2006 representing four seasons. The 

14 derived images for 5 August 2005 and 15 February 2006 were from MODIS-Aqua, while 

15 the other two were from SeaWiFS. 

16 

17 Figure 9. The distribution of the relative percentage of a p h, ad, and a g to their sum at 443 

18 nm within the SMAB. See Figure 8 for detail. 

19 

20 Figure 10. The relative percentage of a p h, ad, and a g to their sum at 443 nm from field 

21 measurements grouped into two regions (near shore region with bottom depth <20 m, and 

22 offshore region with bottom depth >20 m) within the SMAB. 

23 
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1 Figure 1 1 . The distribution of [Chla] calculated from operational ocean color algorithms 

2 (OC4V4 for SeaWiFS and OC3M for MODIS-Aqua) and from the empirical method 

3 described in this paper (OCSMAB; [Chi a] = 70.632 x [ 0 p(I (670)] 1 184 ), and their ratio 

4 [(OC4V4 or OC3M)/OC_SMAB] within the SMAB. See Figure 8 for detail. The 

5 responding scales of the color bar are in log units for [Chl a] and in linear units for the 

6 ratio. 

7 

8 Figure 12. Monthly time series of a) [Chl a], b) a g ( 443), and c) diffuse attenuation 

9 coefficient at 490 nm (K490) from MODIS-Aqua Level-3 images (4x4 km resolution) for 

10 a near shore location (75.90W, 36.93N; solid circle) and a middle shelf location (75.30W, 

1 1 36.93N; open circle). [Chl a] and a g (443) are calculated from algorithms developed in 

12 this paper, while K490 is a direct product from the Level-3 images. Monthly river 

13 discharge rates at the mouth of Chesapeake Bay (open triangle; Data sources: 

14 http://waterdata.usgs.gov/nwis/; written communication from Gary Fisher, U.S. 

15 Geological Survey, 17 July 2007) are also shown for comparison. 
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Table 1 . Regression results of phytoplankton absorption coefficient ( a p h ) to a^/,(670) from 
Equation (3): a ph (A) = B 0 (A)[a ph (670)] BlU) . Log-transformation was applied to the data 

and Model II linear regression was adopted to calculate log[Bo(A.)] and Bi(X). The 
selected wavelengths for a p h analysis were the visible bands for SeaWiFS and MODIS- 
Aqua. 


X (nm) 

May-October 

(N=196) 

Bo Bj r 2 

November- April 
(N=51) 

Bo Bj r- 

412 

1.296 

0.835 

0.986 

1.131 

0.849 

0.987 

443 

1.525 

0.843 

0.989 

1.290 

0.848 

0.979 

488 

1.023 

0.846 

0.983 

0.806 

0.821 

0.961 

490 

1.015 

0.851 

0.983 

0.800 

0.825 

0.963 

510 

0.842 

0.911 

0.983 

0.637 

0.856 

0.975 

531 

0.694 

0.983 

0.970 

0.489 

0.875 

0.974 

551 

0.603 

1.047 

0.945 

0.378 

0.893 

0.956 

555 

0.587 

1.067 

0.938 

0.347 

0.895 

0.948 

667 

0.899 

1.005 

1.000 

0.923 

1.012 

1.000 

678 

1.039 

1.002 

0.999 

0.914 

0.959 

0.994 
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Table 2. Statistical results for absorption coefficients of phytoplankton ( a p h ), non- 
pigmented particles (a d ), and CDOM plus non-pigmented particles (a dg ) at selected 
wavelengths from loglinearmodel: log[a, (/l)] = C 0 (A) + Cj (A)R , where 

R = log[R rs (/lj ) / R rs ( A 2 )] . The size of the dataset is N=25. 

Parameter f? n (490)/f?„(555) f?„(490)/f?„(551) 


Q Co r 2 Cj C 0 r 2 


a^(670) 

-2.602 

-1.467 

0.921 

-2.769 

-1.487 

0.923 

0^380) 

-2.797 

-1.319 

0.933 

-2.976 

-1.340 

0.932 

M 400) 

-2.812 

-1.387 

0.933 

-2.992 

-1.408 

0.931 

a<i 412) 

-2.849 

-1.427 

0.933 

-3.031 

-1.449 

0.931 

ad( 443) 

-3.048 

-1.633 

0.925 

-3.243 

-1.656 

0.924 

0^490) 

-3.260 

-1.950 

0.902 

-3.468 

-1.975 

0.902 

a<( 510) 

-3.584 

-2.105 

0.864 

-3.813 

-2.132 

0.864 

aj( 531) 

-3.676 

-2.223 

0.848 

-3.911 

-2.251 

0.848 

ct<( 555) 

-3.315 

-2.297 

0.891 

-3.526 

-2.322 

0.891 

ad g ( 380) 

-1.394 

-0.434 

0.919 

-1.487 

-0.445 

0.919 

0^(400) 

-1.489 

-0.576 

0.921 

-1.588 

-0.587 

0.920 

^(412) 

-1.535 

-0.651 

0.922 

-1.637 

-0.663 

0.922 

^(443) 

-1.593 

-0.879 

0.923 

-1.698 

-0.891 

0.923 

a dg { 490) 

-1.649 

-1.190 

0.910 

-1.758 

-1.203 

0.910 

o^(510) 

-1.725 

-1.308 

0.899 

-1.840 

-1.322 

0.899 

^g(531) 

-1.633 

-1.406 

0.885 

-1.743 

-1.419 

0.886 

^(555) 

-1.983 

-1.564 

0.822 

-2.117 

-1.580 

0.823 
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Table 3. Statistical results for CDOM absorption coefficient ( a g ) at selected wavelengths 


from the reverse exponential model: 


= H o W + (V ex P [-^2 VK W] • The 

rs \^2 ) 


size of the dataset is N=34. 


Parameter 

H 0 

f?„(490)//? n (555) 

Hi H 2 r 2 

H 0 

fl„(490)/fl„(551) 
Hi H 2 

r 2 

‘ 355) 

0.538 

3.149 

3.978 

0.882 

0.546 

2.805 

3.844 

0.879 

«g(380) 

0.534 

3.015 

6.110 

0.902 

0.542 

2.692 

5.909 

0.900 

a g (400) 

0.540 

2.940 

8.656 

0.915 

0.547 

2.625 

8.366 

0.912 

^412) 

0.523 

2.849 

9.914 

0.916 

0.531 

2.551 

9.592 

0.914 

ag (443) 

0.531 

2.857 

17.700 

0.908 

0.539 

2.557 

17.130 

0.906 

a g { 490) 

0.547 

3.138 

39.960 

0.891 

0.555 

2.798 

38.690 

0.890 

^510) 

0.493 

2.352 

39.870 

0.857 

0.503 

2.126 

38.640 

0.858 

**(531) 

0.494 

2.271 

50.240 

0.829 

0.504 

2.056 

48.700 

0.831 

**(555) 

0.335 

1.798 

40.690 

0.795 

0.346 

1.657 

39.120 

0.799 
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Table 4. The mean absolute percent difference (MAPD) and root mean square error 
(RMSE) from validation match-ups for a p h, dj, a g , and aj g at selected wavelengths. The 
derived exponential decay slope ( S) for ad (Sj), a g (S g ), and ad g (Sd g ) from non-linear 
regression and chlorophyll a concentration ([Chla]) are also compared. Data used for 
algorithm development are not included in this analysis. The size of the datasets are 
N=22, 36, 8, and 19 for a p h, ad, or ad g , and N=29, 45, 14, and 25 for [Chl a], and N=31, 
47, 14, and 25 for a g for SeaWiFS ±8 and ±32 hours and MODIS ±8 and ±32 hours 
overpass windows, respectively. 

Parameter SeaWiFS (±8h) SeaWiFS (±32h) MODIS (±8h) MODIS (±32 h) 


MAPD RMSE MAPD RMSE MAPD RMSE MAPD RMSE 


a ph{ 443) 

23.5 

0.0501 

23.4 

0.0470 

23.4 

0.0137 

23.1 

0.0206 

a ph (6 70) 

25.7 

0.0371 

30.1 

0.0329 

21.2 

0.0051 

27.8 

0.0106 

[Chla] 

32.3 

3.6808 

32.5 

3.0339 

28.8 

3.2805 

29.3 

2.5474 

cid( 380) 

35.3 

0.1447 

38.2 

0.1350 

41.9 

0.1684 

39.3 

0.1275 

ad( 443) 

41.8 

0.0787 

42.1 

0.0726 

50.5 

0.0876 

44.6 

0.0664 

s d 

13.4 

0.0018 

11.7 

0.0016 

11.5 

0.0014 

12.1 

0.0015 

a g ( 380) 

25.2 

0.1636 

20.7 

0.1359 

20.1 

0.2096 

23.4 

0.1607 

%(443) 

22.8 

0.0516 

20.1 

0.0432 

20.1 

0.0662 

21.5 

0.0509 

Sg 

5.8 

0.0011 

5.5 

0.0011 

5.2 

0.0011 

6.7 

0.0012 

ctd g { 380) 

26.5 

0.2104 

25.2 

0.1965 

18.4 

0.1761 

25.6 

0.1478 

a^443) 

24.4 

0.0923 

24.4 

0.0870 

22.1 

0.0896 

23.1 

0.0711 

Sdg 

11.9 

0.0020 

12.7 

0.0021 

15.5 

0.0024 

15.1 

0.0023 
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Longitude 

Figure 1. Map of the study area within the southern Middle Atlantic Bight (SMAB). 
Symbols representing the sampling stations from the following cruises are: A - BIOME1 
(30 March to 1 April 2005), V - BIOME2 (26 to 30 July 2005), O - BIOME3 (5 to 9 
May 2006), X - BIOME4 (2 to 6 July 2006), □ - CBP (four daily cruises), and ■ - 
CBFI (ten daily cruises). 
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Figure 2. Log-transformed linear regression of phytoplankton absorption coefficient at 
670 nm [^(670)] to a) a p h(443), and b) chlorophyll a concentration [Chi a]. The solid 
and dashed lines in the upper figure represent the regression for summer-fall season (May 
to October) and winter-spring season (November to April), respectively. 
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Figure 3. Absorption algorithms derived from field observations of remote sensing 


reflectance ( R rs ) from a) log linear model for a p i,(6 70): log [a pll (X)] = C 0 (X) + C, (X)R , 


where R = log[i? rs (^1, ) / R rs (X 2 )] ; and b) log_polynomial model for non-pigmented 


particulate absorption coefficient at 380 nm [<7^(380)]: 

log[u rf {X)] = D 0 ( X ) + Dj (X)R + D 2 (X )R 2 + D 3 (X)R 3 + D 4 (X)R 4 ; and c) 

exponential model for CDOM absorption coefficient at 380 nm [a g (380)]: 

a g ( X ) = G 0 (X) + Gj (A)exp[-G 2 (X) — — ^ ] ; and d) reverse exponential model for 

Ks^i) 

absorption coefficient by non-pigmented particles plus CDOM at 380 nm [<7^(380)] 

R (A ) 

— - — — = H 0 ( X ) + H x (X) exp [~H 2 ( X)a dg ( X )] . Regression lines from R rs band ratios of 

KMi) 
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ratio log a d (380) 





412/555, 443/555, and 490/555 are represented as solid, dotted, and dashed lines, 
respectively. 
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Figure 4. Validation results comparing SeaWiFS observations with field measurements of 
a p h(61Q), ad( 380), a g (380), and 380) from multiple models (log linear model, 

log_polynomial model, exponential model, reverse exponential model). Figures a) and 
b) show the mean absolute percent difference (MAPD), while c) and d) show the root 
mean square error (RMSE) of the validation results within 8 hours and 32 hours of the 
satellite overpass, respectively. The data from stations applied to develop the algorithms 
were not included in this analysis. The satellite derived o^(670) and a f / g (380) from 
GSM01 model [Maritorena et al, 2002] and GSM01-CB model [Magnuson et al, 2004] 
are also shown for comparison. 
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Figure 5. Validation results comparing MODIS-Aqua observations with field 
measurements of a p h(61Qi), a<i( 380), a^(380), and a</ g (380) from multiple models. See 
Figure 4 for details. 
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field measurements field measurements 

Figure 6. Comparisons of SeaWiFS and MODIS-Aqua and field observations of the 
absorption coefficients of a) a p h{ 670), b) ajQS 0), c) a g ( 380), and d) <2^(380). The values 
are plotted on log scale. The satellite derivations of a p h{ 670), ad(380), and ^^380) were 
from the loglinearmodel, while <3^(3 80) from the reverseexponentialmodel. The 
match-ups procedure is limited to within ±8 hours, and the data from stations used to 
develop algorithms are excluded for validation analyses. The statistical results are based 
on log-transformation of the data and shown on upper left for SeaWiFS and lower right 
for MODIS-Aqua. The solid lines represent the 1:1 lines, while dashed lines and dotted 
lines represent the regression for SeaWiFS and MODIS-Aqua respectively. 
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Figure 7. Comparisons of SeaWiFS and MODIS-Aqua and field observations of [Chl a] 
for satellite overpass window of a) ±8 hours and b) ±32 hours. The data from stations 
used to develop algorithms of a p h are excluded for this analysis. The solid lines represent 
the 1 : 1 lines, while dashed lines and dotted lines represent the regression for SeaWiFS 
and MODIS-Aqua respectively. See Figure 6 for detail. 
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Figure 8. The distribution of a p h, ad, and a g at 443 run within the SMAB for 5 August and 
3 November 2005, and 15 February and 12 May 2006 representing four seasons. The 
derived images for 5 August 2005 and 15 February 2006 were from MODIS-Aqua, while 
the other two were from SeaWiFS. 
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Figure 9. The distribution of the relative percentage of a p h, a d, and a g to their sum at 443 
nm within the SMAB. See Figure 8 for detail. 
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Figure 10. The relative percentage of a p h, and a g to their sum at 443 nm from field 
measurements grouped into two regions (near shore region with bottom depth <20 m, and 
offshore region with bottom depth >20 m) within the SMAB. 
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Figure 1 1 . The distribution of [Chla] calculated from operational ocean color algorithms 
(OC4V4 for SeaWiFS and OC3M for MODIS-Aqua) and from the empirical method 
described in this paper (OC SMAB; [Chi _a] = 70.632 x [a plt (670)] 1 184 ), and their ratio 


[(OC4V4 or OC3M)/OC_SMAB] within the SMAB. See Figure 8 for detail. The 
responding scales of the color bar are in log units for [Chl a] and in linear units for the 
ratio. 
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Monthly time series 

Figure 12. Monthly time series of a) [Chla], b) 443), and c) diffuse attenuation 
coefficient at 490 nm (K490) from MODIS-Aqua Level-3 images (4x4 km resolution) for 
a near shore location (75.90W, 36.93N; solid circle) and a middle shelf location (75.30W, 
36.93N; open circle). [Chl a] and a g { 443) are calculated from algorithms developed in 
this paper, while K490 is a direct product from the Level-3 images. Monthly river 
discharge rates at the mouth of Chesapeake Bay (open triangle; Data sources: 
http://waterdata.usgs.gov/nwis/; written communication from Gary Fisher, U.S. 
Geological Survey, 1 7 July 2007) are also shown for comparison. 
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